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Abstract 

The local search framework for obtaining PTASs for NP-hard geometric optimiza- 
tion problems was introduced, independently, by Chan and Har-Peled [Hj and Mustafa 
and Ray |17j . In this paper, we generalize the framework by extending its analysis to 
additional families of graphs, beyond the family of planar graphs. We then present 
several applications of the generalized framework, some of which are very different 
from those presented to date (using the original framework). These applications in- 
clude PTASs for finding a maximum l-shallow set of a set of fat objects, for finding a 
maximum triangle matching in an ^-shallow unit disk graph, and for vertex- guarding a 
(not-necessarily-simple) polygon under an appropriate shallowness assumption. 

We also present a PTAS (using the original framework) for the important problem 
where one has to find a minimum-cardinality subset of a given set of disks (of varying 
radii) that covers a given set of points, and apply it to a class cover problem (studied 
in [3]) to obtain an improved solution. 

1 Introduction 

In their break-through papers, Chan and Har-Peled [6] and Mustafa and Ray [T7] showed, 
independently, how a simple local-search-based algorithm can be employed to obtain a 
PTAS for an NP-hard geometric optimization problem. Chan and Har-Peled [6] used 
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local search to obtain a PTAS for finding a maximum independent set of pseudo disks, 
and Mustafa and Ray [T7] used it to obtain a PTAS for finding a minimum hitting set 
(from a given set of points) for half-spaces in M 3 and for r-admissible regions in M 2 . This 
technique turned out to be very powerful, and since its publication, it was applied to a 
variety of additional problems. Gibson et al. [13] used it to obtain a PTAS for the 1.5D 
terrain guarding problem, and Gibson and Pirwani used it to obtain a PTAS for finding a 
dominating set in disk graphs [14] . 

The local-search-based algorithm, as described in [6|I17]. receives an integer parameter 
k and proceeds as follows. It starts with any feasible solution and performs a series of local 
improvements, where each such improvement involves only 0(k) objects. The analysis 
relies on the existence of a planar bipartite graph G, whose vertices on one side correspond 
to the objects found by the local search algorithm ("blue vertices") and on the other side 
to the objects in an optimal solution ("red vertices"), and whose ("blue-red") edges satisfy 
an appropriate locality property relating the two solutions. Chan & Har-Peled and Mustafa 
& Ray showed that for the problems mentioned above such a planar graph G exists, and 
applied the planar separator theorem to relate the size of the local search solution with 
that of the optimal solution. 

In this paper we generalize the local search technique by extending its analysis to 
additional families of graphs. We show that one can achieve a PTAS even in cases where 
the graph whose vertices are the elements of the two solutions and which satisfies the 
locality property, is not planar, but belongs to a family of graphs that has a separator 
property. It is well known that there are many such families of graphs, e.g., graphs with 
forbidden minors [1] and various intersection graphs [H)l[l6j[18]. We also present several 
interesting applications of our extended analysis, some of which are very different from 
those presented in the past (using the original analysis). 

These applications, which are presented in Section O include finding a maximum l- 
shallow set of a set of fat objects, finding a maximum triangle matching in a unit disk 
graph, and guarding a polygon with limited range of sight. We briefly discuss each of them. 
Maximum l-shallow set. Let D be a set of n fat objects in ]R 2 and let I > be a constant. 
An l-shallow subset of D is a subset S of D, such that the depth of the arrangement of S is 
at most /, i.e. every point in the plane is covered by at most I objects of S. In the maximum 
/-shallow subset problem, one is asked to find a maximum-cardinality /-shallow subset S of 
D. Notice that for I = 1, S is a maximum independent subset of D. Chan [1] presented a 
PTAS for the problem of finding a maximum independent subset of D. Later, as mentioned 
above, Chan and Har-Peled [6] presented a local-search-based PTAS for finding a maximum 
independent set of pseudo-disks (and also of D). Notice that a maximum /-shallow set can 
be larger than the set that is obtained by repeatedly finding a maximum independent set 
of the set of remaining objects. We show that our generalized analysis enables us to obtain 
a PTAS for the maximum /-shallow subset problem, and emphasize that it is essential 
whenever / > 1. 

Maximum triangle matching. Given a graph G, one has to find a maximum-cardinality 
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collection of vertex-disjoint triangles in G. Baker [2] presented a PTAS for the important 
case where G is planar. We give a PTAS for the case where G is an /-shallow unit disk 
graph. This PTAS also holds in more general settings, see below. 

Guarding. We apply our generalized analysis to several guarding problems. Many guard- 
ing problems are known to be APX-hard (e.g., guarding a simple n-gon P with as few 
vertex guards as possible [9]), and as such do not admit a PTAS. Ghosh [12] presented 
an 0(log n)-approximation algorithm for vertex-guarding a polygon. Subsequently, Efrat 
and Har-Peled [8] presented an 0(log |OPT|)-approximation algorithm for this problem. 
In contrast to these results, we obtain a PTAS for this problem, under the assumption that 
each vertex guard g has a limited range of sight r g , such that every point in P is covered 
by at least one of the guards and the set of disks centered at the guards is /-shallow. We 
also discuss several other versions, including guarding through walls and guarding a 1.5D 
terrain with bounded range of sight. 

Finally, in Section HI we consider the important problem known as discrete coverage of 
points by disks, for which we obtain a PTAS within the original local search framework. 
Let P be a set of points in the plane and let D be a set of disks that covers P. One needs 
to find a minimum-cardinality subset D' C D that covers P. Notice that the special case 
where D is a set of unit disks, is the dual of the discrete hitting set problem for unit disks, 
for which Mustafa and Ray [17] presented a PTAS (actually, for arbitrary disks). However, 
in the general case, where D is a set of disks of varying radii, coverage and hitting are not 
dual. We present a PTAS for the general case that is inspired by the work of Gibson and 
Pirwani p3j. We also apply this result to the a class cover problem, improving a result of 
Bereg et al. [3]. 

2 PTAS via local search 

We begin by generalizing the local search technique to additional families of graphs, beyond 
the family of planar graphs, such as intersection graphs and graphs with forbidden minors. 
Actually, we describe the technique for any family of graphs that has a separator property, 
similar to the separator property of planar graphs. 

Let J 7 be a monotone family of graphs, i.e., all subgraphs of a graph G € J- are also 
in T. Assume that T has a separator property, i.e., for any graph G = (V, E) in J 7 , one 
can partition the vertex set V, where \V\ = n, into three sets A,B and S, such that (i) 
\S\ < en 1 "" 5 , where < 5 < 1 and c > is a constant, (ii) \A\, \B\ < an, where 1/2 < a < 1, 
and (iii) the sets A and B are disconnected, i.e., there is no edge in E between a vertex of 
A and a vertex of B. 

Frederickson |11] defined the notion of an r-division for planar graphs, and showed 
how to obtain an r-division for a given planar graph G, by repeatedly applying the pla- 
nar separator algorithm. It is straight-forward to adapt Frederickson's construction and 
analysis to the family T . Essentially, one needs only to replace the exponent 1/2 in the 
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size of the separator by 1 — 5. Thus, J- has the following property. Let r be a parameter, 
1 < r < n, then for any connected graph G = (V, E) in T , one can find a collection of 
Q(n/r) pairwise disjoint subsets V±, V2, ■ ■ ■ of V, such that (i) \V,\ < C2r, where C2 > is a 
constant, (ii) |r(Vj)| < c^r 1-5 , where T(Vi) is the set of vertices in V \ V{ that are adjacent 
to a vertex in V, and C3 > is a constant, and (iii) UT(V) = S, where S = V \ UV^; in 
particular, T(Vi) C S, and, for any j ^ i, Vi and Vj are disconnected. It follows that 
IS"! < Q(n/r) ■ C3r 1-5 = c\^, where c\ > is a constant. Also, it is easy to verify that 

Claim 2.1. If r = C2 ^_ C3 , where k is a sufficiently large constant, then, for any index i, 

W\ + |r(^)| < k. 

We distinguish between minimization problems and maximization problems. 
2.1 Minimization problems 

Let us now recall the local search technique as it is used, e.g., in the context of geometric 
piercing or covering. Assume that we are considering a minimization problem V, that is, 
one needs to find a minimum-cardinality subset of a given set X that is a solution for V. 
Let Xq C X be an initial not-necessarily-optimal solution for V, and let k be a sufficiently 
large constant. The local search technique checks whether there exists a subset X' C Xq 
of size k and a subset X" C X of size k — 1, such that (Xq \ X') U X" is still a solution for 
"P. If yes, then it replaces X' by X" (i.e., it performs this local improvement) and resumes 
the search. Otherwise, it halts. 

Let B and R be two solutions for V, where B was obtained by applying local search, 
and assume that B n R = 0. (Otherwise, we can remove the elements that belong to 
both B and R.) Moreover, let G = (V,E) be a graph, such that (i) V = B U R, and (ii) 
for each object o that needs to be "solved", there exists an edge e £ E between a vertex 
b € B and a vertex r £ R, where both b and r "solve" o. This requirement is called the 
locality condition. We prove below that if G belongs to some family T that has a separator 
property, then |J5| is not much greater than \R\. Parts of the proof appear already in [6lll7j 
and are included here for completeness. We first construct an r-division of G, for r = C2 ^_ C3 , 
and set Bj = B n Vi and R4 = Rn Vi. 

Claim 2.2. For any index i, the set (B \ Bi) U T(Bi) is also a solution for V. 

Proof. Fix i and let o be an object that needs to be "solved". If all vertices of B that 
"solve" o belong to Bi, then, by the locality condition, there exists r € T(Bi) that "solves" 
o. Otherwise, there is a vertex b € B \ Bi that "solves" o. We conclude that in both cases 
there is a vertex in (B \ Bi) U V(Bj) that "solves" o. □ 

Claim 2.3. For any index i, \Bi\ < \Ri\ + \T(V)\ < k. 
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Proof. Observe first that by arguments similar to those in the proof of Claim 12.21 the set 
(B \ Bi) U {Ri U T(Yi)) is also a solution for V . Moreover, by Claim EH \Ri\ + |T(T^)| < 
\Vi\ + |r(Vj)| < k. So, if \Bj\ > \Ri\ + |r(Vi)|, then the local search algorithm would have 
replaced Bi (or a subset of Bi of size k, if \Bi\ > k) by R4 U r(T^) before halting. Since it 
has not done so, we conclude that \Bi\ < \Ri \ + |T(Vi)|. □ 

Theorem 2.4. For any e, < e < 1, one can choose a constant k, such that \B\ < 
(l + e)\R\. 

Proof. Set d = 2ci(c 2 + c 3 ) 5 , e' = §, k = and recall that r = ^J^. Then, \B\ < 

|S|+E, \Bi\ < l^l+E, l^|+Ei \nVi)\ < ™ + \R\ + w = \R\ + *zp = \R\+\*£W2ci(c 2 + 
c 3 ) 5 = \R\+e'(\R\ + \B\). We thus have, \B\ < \R\ ■ i±|^ = (1 + e')(l + e' + e' 2 + . . .) < 
\R\(1 + e')(l + e' + e') = \R\(1 + 3e' + 2e' 2 ) < |i?|(l + 5e') = |fl|(l + e). 

□ 

2.2 Maximization problems 

Let us now recall the local search technique as it is used, e.g., in the context of geometric 
packing. Assume that we are considering a maximization problem V, that is, one needs to 
find a maximum-cardinality subset of a given set X that is a solution for V . Let Xq C X be 
an initial not-necessarily-optimal solution for V, and let k be a sufficiently large constant. 
The local search technique checks whether there exists a subset X' C of size at most 
fe — 1 and a subset X" C X of size k, such that (Xo \ X') U X" is still a solution for P. If 
yes, then it replaces X' by X" (i.e., it performs this local improvement) and resumes the 
search. Otherwise, it halts. 

Let B and R be two solutions for V, where B was obtained by applying local search, 
and assume that B n R = 0. (Otherwise, we can remove the elements that belong to both 
B and R.) Moreover, let G = (V, E) be a graph, such that (i) V = B L) R, and (ii) there 
exists an edge e € E between a vertex b € -B and a vertex r S i? if and only if b and r 
intersect. 

Theorem 12.71 below states that if G belongs to some family T that has a separator 
property, then \B\ > (1 — e)|i?|. The proof is similar to the one in Section r2.lt we include it 
for completeness. We first construct an r-division of G, for r = c _^ C3 , and set Bi = B n ^ 
and i?i = i?n Vi. 

Claim 2.5. For any index i, the set (B \ (Bi U T(Vi))) U Ri is also a solution for V . 

Proof. Fix i, and let r £ Ri. By definition, r intersects only its neighbors in G, i.e., only 
elements in (Bi U r(V5)) U R4. Thus, for each b E B \ (Bi U r(T^))), r n 6 = 0. Therefore, 
(5 \ (Bi U r(Vi))) U Ri is also a solution for V. □ 

Claim 2.6. For any index i, \Ri\ < \Bi\ + \T{y.j)\ < k. 
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Proof. By Claim [231 the set (* \ (Bi U T(Vi))) U Ri is also a solution for "P. Moreover, by 
ClaimEU |B<ur(Vi)| = |J3i| + |r(^)| < |^| + |r(V-)| < ft. So, if \Ri\ > \BiUT(Vi)\, then the 
local search algorithm would have replaced Bi U T(Vi) by Ri (or a subset of Ri of size k, if 
|*| > k) before halting. Since it has not done so, we conclude that \Ri\ < \Bi\ + \T(Vi)\. □ 

Theorem 2.7. For any e, < e < 1, one can choose a constant k, such that \B\ > 
(l-e)|*|. 

Proof. Set d = 2ci(c 2 + C3)* 5 , k = (c'(f - l)) 1 / 5 , and recall that r = Then, 

|*| < |5| + £ 1*1 < \S\ + E 1^1 + E I W)l < ^ + \B\ + ^ = |B| + ^ 
= 1^1 + |jR| ^ |jB| 2ci(c 2 + c 3 ) 5 = |*| + j^(\R\ + \B\) = \B\ + 237(1*1 + 1*1) • 

e 

Rearranging, we get that |*| > (1 — e)|*|. □ 

3 Applications 

In this section we describe several original applications of our generalized local search 
technique. The applications in Sections 13.11 and 13.21 use the maximization version and the 
applications in Section [3.31 use the minimization version. 

3.1 Maximum /-shallow set for fat objects 

We consider the maximum I -shallow set problem for a set of fat objects. Let * be a set 
of n fat objects and let I > be a constant. An I -shallow set of * is a subset 5 of *, 
such that the depth of the arrangement of S is at most I. In the maximum /-shallow set 
problem, one has to find a maximum-cardinality subset S of *, such that S is /-shallow. 
Notice that for / = 1, 5 is a maximum independent set of *. 

We show that our generalization of the local search technique enables us to apply local 
search to find a (1 — e)-approximation of a maximum /-shallow set of *. Indeed, let * be 
the set of fat objects obtained by applying local search and let * be an optimal set. (We 
may assume that * n * = 0, since otherwise, we can remove the objects that appear in 
both sets.) Notice that the analysis of Chan and Har-Peled [6] does not immediately apply 
here (assuming / > 1), since the intersection graph of *U* is not necessarily planar (even 
if one removes the "monochromatic" edges) However, this graph does have a separator of 
size 0{^/n). This follows from results of Miller et al. [16] and Smith and Wormald |18j . 
who show that the intersection graph of a set of /-shallow fat objects has a separator of 
size v/n. 
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Recall that the local search algorithm begins with the empty solution, and in each 
iteration it improves the current solution by replacing a subset of size at most k — 1 of 
the current solution with a larger subset of size at most k of D, such that the resulting 
set is still a solution (i.e., it is still /-shallow). Since k is a constant, each iteration can 
be performed in polynomial time, and therefore the total running time of the local search 
algorithm is polynomial in n. 

As for the size of B, since both B and R are /-shallow, the set B U R is 2/-shallow, and 
therefore the intersection graph of B U R has a separator of size \/lln = 0(^/n). Hence, 
by Theorem 12.71 we conclude that \B\ > (1 — e)|-R|. 

Finally, we note that Chan's separator-based method [4] can be used to obtain a (1 — e)- 
approximation of a maximum /-shallow set of D. However, as also mentioned in [5J, the 
disadvantage of his method is that it explicitly applies an algorithm for finding a separator, 
while we apply it only for analysis purposes. 

3.2 Maximum triangle matching in unit disk graphs 

Next, we consider the maximum triangle matching problem in unit disk graphs. Given 
a graph G = (V,E), find a maximum-cardinality collection of pairwise-disjoint subsets 
Vi, . . . , V m of V, each consisting of exactly 3 vertices, such that for each Vi = {ui,Vi, Wi}, 
1 < i < m, the three edges (ui, v^), (yi, Wi), (wi, Ui) belong to E. Baker [2] presented a 
PTAS for the important case where G is planar. We give a PTAS for the case where G is 
an /-shallow unit disk graph. This PTAS also holds in more general settings; see remark 
below. 

Let S be a set of n points in the plane, and let UDG(S) be the graph over S, in which 
there is an edge between u and v if and only if the Euclidean distance between u and v is 
at most 1. Let V be the set of disks of radius 1/2 centered at the points of S, and assume 
that the depth of the arrangement of T> is at most some constant /. Then, there is an edge 
between u and v if and only if d u n d v ^ 0, where d u , d v are the disks of V centered at u, v, 
respectively. 

Our generalization of the local search technique enables us to apply local search to find 
a (1 — e)-approximation of a maximum triangle matching of UDG(S). We begin with the 
empty set of triangles, and in each iteration we replace a subset of at most k — 1 triangles 
of the current solution with a larger subset of at most k triangles, such that the resulting 
set is still a solution (i.e., it is still a triangle matching). 

Consider the set of triangles B obtained by applying local search, and the set of triangles 
1Z of a maximum matching. (Each A £ B U 1Z is a triangle A = {a, b, c}, where a, b, c are 
points in S and (a, 6), (b, c), (c, a) are edges of UDG(S).) For each triangle A in B (resp. 
in TV), select an arbitrary unique representative point p/\ that lies in the interior of A, and 
denote by B (resp. R) the resulting set of representative points. 

Now, consider the graph G = (B U R, E), where (p, q) £ E if and only if the Euclidean 
distance between p and q is at most 2. Notice that G satisfies the following locality 
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condition. For any two points b £ B and r 6 fi, whose corresponding triangles are 
AubVbWt and Au r v r w r , respectively, and such that {ub, Vb, Wb} D {ii r , « r , tt! r } 7^ 0, the edge 
(b, r) belongs to E. This is true since the distance from b (alt., r) to any of the its triangle 
corners is at most 1. 

It remains to show that G has a separator. Let T>' be the set of unit disks centered 
at the points of B U R. Then, (p, q) £ E if and only if the disks of V centered at p and 
q, respectively, intersect. It is easy to see that the depth of the arrangement of T>' is at 
most some constant c = c(l). This follows from the assumption that the depth of the 
arrangement of T> is at most I. Therefore, by Miller et al. |16j . G has a separator of size 
0(y/cn) = 0(y/n). Hence, by Theorem 12.71 we conclude that \B\ > (1 — e)|i?|. 
Remark. It is easy to see that one can replace "triangle" in the above result with any 
connected graph (i.e., structure) H with a constant number of vertices, and obtain a PTAS 
for maximum //-matching in unit disk graphs. 

3.3 Guarding with limited visibility 

We now demonstrate our generalization of the local search technique on a family of covering 
problems, or, more precisely, on a family of guarding problems. Throughout this section, 
Q denotes a set of points, representing stationary guards, where each guard g £ Q has 
its own range of sight r g . Let D g denote the disk of radius r g centered at g, and set 
D = {D g : g £ Q}. We assume that the depth of the arrangement of D is at most some 
constant I. 

Guarding a polygon. Given a polygon P (possibly with holes) and a set Q C P, a 

minimum guarding set for P (with respect to Q) is a minimum-cardinality subset Q' of Q, 
such that every p £ P is guarded by Q' , i.e., for every p £ P, there exists g £ Q' , such that 
\gp\ < r g and the segment gp is contained in P. 

Eidenbenz et al. [9] proved that finding a minimum guarding set for a polygon, where 
the given set of guards is the polygon's set of vertices, is APX-hard, even if the polygon 
has no holes and there are no limits on the ranges. We show that our generalization of the 
local search technique enables us to apply local search to find a (1 + e)-approximation of a 
minimum guarding set for P (with respect to Q), under the above assumption concerning 
the depth of the arrangement of D. 

For a guard g £ Q, let VP g denote its visibility polygon and let VR g denote its visible 
region, where VR g = VP g n D g . Notice that VR g is not necessarily convex. For a subset 
X C Q, let D x denote the set {D g : g £ X} and let VR X denote the set {VR g : g £ X}. 
Let B C Q be the set of guards obtained by applying local search (to the set VRg) and let 
R be an optimal set. (We may assume that B n R = 0, since otherwise, we can remove the 
guards that appear in both sets.) Consider the graph G = {B U R,E), where (g,g') £ E 
if and only if VR g D VR g > / 0. We claim that G has a separator of size \fhi = 0(y/n), 
since any separator of the intersection graph of DgU Dr is also a separator of G, and the 
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former graph has a separator of size Vm = 0(y/n), by a result of Miller et al. [16]. Observe 
that the locality condition is satisfied, since for any p £ P, let b £ Z? and r £ i? be two 
guards that guard p, then VRft n V7? r 7^ and hence the edge (b, r) is in G. Therefore, by 
Theorem 12. 4| we conclude that < (1 + e)|i?|. 

Theorem 3.1. There exists a PTAS for minimum guarding a (not necessarily simple) 
polygon with respect to a given set of guards, assuming that the depth of the corresponding 
arrangement of disks is bounded by a constant. 

It is possible to consider other versions of the problem in which the visible region of a 
guard g £ Q is defined differently, as long as VR g C D g . The approximation analysis for 
these versions is the same as for the standard version. For example, consider the problem 
of guarding a polygon P, where each guard can see through at most some fixed number of 
walls (edges of the polygon P). Then, the visibility polygon of a guard g £ Q is different 
than its visibility polygon in the standard version, but still its visible region is contained in 
D g . For such problems, one needs to modify the verification procedure of the local search 
algorithm (i.e., whether a given set of guards consists of a solution). We thus obtain, for 
example, the following corollary. 

Corollary 3.2. There exists a PTAS for minimum guarding a polygon through walls, 
assuming that the depth of the corresponding arrangement of disks is bounded by a constant. 

Terrain guarding with limited visibility. Let T be a 1.5D terrain (i.e., an x- monotone 
polygonal chain), let Q C T be a finite set of guards, and let X C T be a finite set of points 
to be guarded. A guard g £ Q sees a point x £ X if ~gx does not cross any edge of T and 
\gx\ < Tg. The goal is to find a minimum-cardinality subset Q' C Q, such that Q' sees X 
(i.e., for each point x £ X, there exists a guard in Q' that sees x). Assuming unlimited 
visibility (i.e., r g = 00, for each g £ Q), the problem is known to be NP-hard [HQS], an d 
there exists a local-search-based PTAS for it [13]. This PTAS relies heavily on the, so 
called, "order claim" , which states that for any four points a, b, c, d on T (listed from left 
to right), if a sees c and b sees d then a also sees d. Unfortunately, this claim is false in 
the limited visibility version. However, this version is a simple variant of polygon guarding 
with limited visibility. Therefore, we have 

Theorem 3.3. There exists a PTAS for terrain guarding, assuming that the depth of the 
corresponding arrangement of disks is bounded by a constant. 

Consider the guarding problems above, but now set r g = 1, for each g £ Q, and assume 
the requirement is to guard a finite set X of points in P (alt., on T). Then, one may replace 
the assumption that the depth of the arrangement of disks is bounded by a constant with 
the assumption that the points in X are not too crowded, i.e., the distance between any two 
points in X is at least some constant 5 > 0. This assumption is sometimes more convenient. 
Under this assumption, each guard sees at most some constant number c = c(5) of points of 
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X, and it is likely therefore that there are many redundant guards. (A guard is redundant 
if there exists another guard that sees the same subset of points of X.) We thus remove 
redundant guards from Q, one at a time, as long as there are such guards in Q. It is easy 
to see that the depth of the arrangement of disks corresponding to the set of remaining 
guards is bounded by a constant. We obtain, e.g., the following theorem. 

Theorem 3.4. There exists a PTAS for minimum guarding a sparse set of points X within 
a (not necessarily simple) polygon P with respect to a given set of guards of unit visibility 
range. 

4 Discrete coverage of points 

We consider the following fundamental problem. Let P be a set of points in the plane and 
let D be a set of disks that covers P. (We assume that the centers of the disks in D are in 
general position.) Find a minimum-cardinality subset D' C D that covers P. 

If D is a set of unit disks, then this is the dual version of the problem known as discrete 
piercing of unit disks. Mustafa and Ray [T7] presented a PTAS for the latter problem (even 
for the case of varying radii) . We present a local-search-based PTAS for discrete coverage 
of points. Our proof is based on the framework developed in [T7] and in Section [21 and is 
inspired by the work of Gibson and Pirwani |14j . 

Theorem 4.1. There exists a PTAS for discrete coverage of points by disks (alt., by axis- 
parallel squares). 

Proof. We may assume that there do not exist two disks in D, such that one of them is 
contained in the other, since, if two such disks do exist, we can always prefer the larger 
one. Let B be the set of the centers of the disks returned by the local search algorithm, 
and let R be the set of the centers of the disks in an optimal solution. We may assume 
that B f]R = (since otherwise, we can remove the centers that belong to both sets). For 
a center c £ B U R, let D c £ D denote the corresponding disk and r c its radius. In order 
to prove that \B\ < (1 + e)|-R|, we need to present a planar bipartite graph G = (B U R, E) 
satisfying the following locality condition: For each p £ P there exist b £ B and r £ R, 
such that p £ Db, p £ D r , and (6, r) £ E. 

Consider the additively weighted Voronoi diagram of B U R, constructed according to 
the distance function S(p,c) = d(p,c) — r c . Observe first that for any center c £ B U R 
and point p £ M 2 , if 5(p, c) < 0, then p £ D c (and vice versa). Now, let cell(c) denote the 
cell of the diagram corresponding to c £ B U R. It is well known that for each c £ B U R, 
c £ cell(c) and cell(c) is connected. 

Our planar bipartite graph is the dual graph of the weighted Voronoi diagram defined 
above, without the monochromatic edges. That is, there is an edge between two centers c 
and c' in B U R, if and only if c £ B and d £ R, or vice versa, and their cells are adjacent 
to each other. We denote this graph by G = (R U B, E). It is easy to see that G is planar. 
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Indeed, by the diagram's properties, any edge (c, c') E E can be drawn such that it is 
contained in cell{c) U celiac'), and within each cell cell(c), it is easy to ensure that the 
edges do not cross each other. 

We now show that G satisfies the locality condition. Let p E P, and assume w.l.o.g. 
that p E cell(r), where r E R. Let b E B be the closest center to p according to 5. That is, 
for any b' E B, 5(p, b) < 5(p, b'). Notice that 5(p, b) < 0, since there exists a center b' E B 
whose disk covers p, and therefore, 5(p, b') < 0. We conclude by the observation above 
that p E .CV Consider the cells we visit when walking along the line segment pb from p to 
b. Since p ^ cell(b) and b E cell(b), we must at some point enter cell(b). Let c E B U i? be 
the center for which cell(c) is the last cell that we visit before entering cell(b), and let q 
be the point on pb, which is also on the boundaries of cell(c) and cell(b). 

It remains to show that c E R, implying that (b, c) E £7, and that p E -D c . Using 
the triangle inequality and since the centers of the disks are in general position, we get 
that 5(p, c) = d(p, c) — r c < d(p, q) + d(q, c) —r c = d(p, q) + 5(q, c). But, d(p, q) + 5(q, c) = 
d(p,<l) + S(<l,b) = d(p,q) + d(q,b)-r h = d(p,b)-r b = 5(p,b), so we get that 5(p, c) < 5(p,b). 
Now, since b is the closest center to p among the centers in B, we conclude that c E R, 
and since S(p,b) < 0, we conclude that p E D c . The axis-parallel version is obtained by 
replacing the L2 metric by the metric. □ 

4.1 Discrete coverage of a polygon 

Consider now the following problem. Let Q be a polygon, and let D be a set of disks 
(alt., a set of axis-parallel rectangles). Find a minimum-cardinality subset of D that covers 
Q. The local search algorithm can be easily adapted to this setting. At each iteration 
of the algorithm, instead of checking whether all points are covered, one needs to check 
whether the entire polygon is covered. This can be done in polynomial time. The analysis 
is essentially the same as for discrete coverage of points, except that here p is any point in 
the polygon Q. We thus conclude that 

Theorem 4.2. There exists a PTAS for discrete coverage of a polygon by disks (alt., by 
axis-parallel squares). 

4.2 The class cover problem 

The class cover problem is defined as follows: Let B be a set of blue points and let R be a 
set of red points and set n = \B\ + \R\. Find a minimum-cardinality set D of disks (alt., of 
axis-parallel squares) that covers the blue points, but does not cover any of the red points. 
That is, find a minimum-cardinality set D, such that B n Ud^Dd = B and R(~) UdeDd = 0. 
Bereg et al. [3] study several versions of the class cover problem with boxes. In particular, 
they prove that the class cover problem with axis-parallel squares is NP-hard and give 
an 0(l)-approximation algorithm for this version. We show that the class cover problem 
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with squares (resp., disks) is essentially equivalent to discrete coverage of points by squares 
(resp., disks), and therefore there exists a PTAS for both versions. 

Indeed, consider any set D of disks that cover the points in B and does not cover any 
point in R. It is easy to see that one can replace each disk d G D with a "legal" disk d! , 
such that BflciC B n d! and either there are three points on d"s boundary, or there are 
exactly two points on d"s boundary and the line segment between them is a diameter of 
d! . Therefore, we can transform the class cover problem with disks to discrete coverage 
of points by disks, by first computing all 0(n 3 ) disks defined by either triplets or pairs of 
points and then removing those that are "illegal". Similarly, we can transform the class 
cover problem with squares to discrete coverage of points by squares. We thus obtain 

Theorem 4.3. There exists a PTAS for the class cover problem with disks (alt., with 
axis-parallel squares). 

Concluding remark. Very recently it has been brought to our attention that Chan 
and Grant [5] observe that the PTAS of Mustafa and Ray [17] for discrete hitting set 
of half-spaces in R 3 implies a PTAS for discrete coverage of points by disks, by a lifting 
transformation that maps disks to lower half-spaces in R 3 and by duality between points 
and half-spaces. However, our PTAS above for discrete coverage of points by disks is direct 
and refrains from moving to R 3 . 
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